function beaming=getBeaming(model,m_tr)


switch model
    case 1
        
        beaming=...
            73./m_tr.^2.*(m_tr>8.5)...
            +1.*(m_tr<=8.5);
    case 2
        beaming=...
            (1-(1+2.*73./m_tr.^2).^-0.5).*(m_tr>8.5)...
            + 1.* (m_tr<=8.5);
end
end